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Abstract. - In an effort to understand "rafts" in biological membranes, we propose phe- 
nomenological models for saturated and unsaturated lipid mixtures, and lipid-cholesterol mix- 
tures. We consider simple couplings between the local composition and internal membrane 
structure, and their influence on transitions between liquid and "gel" membrane phases. As- 
suming that the gel transition temperature of the saturated lipid is shifted by the presence 
of the unsaturated lipid, and that cholesterol acts as an external field on the chain melting 
transition, a variety of phase diagrams are obtained. The phase diagrams for binary mix- 
tures of saturated/unsaturated lipids and lipid/cholesterol are in semi-quantitative agreement 
with the experiments. Our results also apply to regions in the ternary phase diagram of 
lipid/lipid/cholesterol systems. 



Introduction. - Biological membranes typically contain various components such as lipid 
mixtures, sterols, and proteins that are indispensable to cell functions [1]. Rather than being 
uniformly distributed in the membrane, there is growing evidence that some intra-membrane 
components are incorporated in domains arising from lateral lipid segregation in membranes. 
This phenomena has attracted great interest in the context of "rafts" [2], i.e., liquid do- 
mains rich in cholesterol, saturated lipids (typically sphingomyelin lipids), and in some cases 
particular proteins [3]. Moreover, cholesterol-rich domains have been directly observed in 
model membranes composed of lipid mixtures and cholesterol, using advanced fluorescence 
microscopy [4-10]. 

Prior to the notion of "rafts" in biological membranes, the role of cholesterol was inves- 
tigated in binary lipid-cholesterol membranes [11-13], where phase separation was observed 
using NMR and calorimetry [14, 15]. Lipid membranes undergo a freezing or "gel-like" tran- 
sition, in which the hydrocarbon tails order. Addition of cholesterol has several effects. It 
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suppresses the "gel" transition below physiologically relevant temperatures, and can lead 
to coexistence of two liquid phases with very different orientational order. It is now be- 
lieved that model membranes containing two phospholipids (saturated and unsaturated) and 
cholesterol [4-10] exhibit "rafts" which are liquid- ordered (Lo) domains, coexisting with a 
surrounding background in a liquid-disordered (Ld) state [3]. 

From a physical viewpoint, a strategy for understanding the basic structure of "rafts" in 
biological membranes is as follows [16]. First, it is necessary to have a simple model for binary 
mixtures of saturated and unsaturated lipids. Second, a minimal model describing binary lipid- 
cholesterol systems is required in order to understand the effects of cholesterol on membrane 
phase behavior. Finally, these two viewpoints could be combined to fully investigate three- 
component systems. Here we focus on the first two steps, and propose simple phenomenological 
models for lipid-lipid and lipid-cholesterol binary mixtures. Such an approach is quite useful 
since we can predict, at least qualitatively, the complex phase behavior [shown, e.g., in fig.|3a)] 
of the ternary lipid-lipid-cholesterol system from the three binary sub-systems. 

It is generally believed that the gel phase transition is driven by the freezing of lateral 
motion as well as conformational ordering of lipids. However, because these two degrees of 
freedom are strongly coupled, the gel phase transition occurs at a single temperature for 
pure lipid systems. One of our major assumptions is that the membrane state, even for lipid 
mixtures, can be described by one internal degree of freedom, coupled to the lateral phase 
separation. Our main result is phase diagrams for two component systems, which reproduce 
experimental ones, without specifying the detailed microscopic state of the lipids. Moreover, 
we shed light on recent phase diagrams obtained for ternary mixture [8,9]. 

Model for saturated-unsaturated lipids systems. - If different lipids exhibit only a small 
difference in their gel transition temperature, their phase diagram in terms of temperature- 
composition parameters will include a "cigar-like" shape of liquid-gel coexistence [as in fig.^a)]. 
As the gel transition temperature difference between the two lipids is increased, a gel-gel co- 
existence region appears below the liquid-gel coexistence. At even larger difference, the phase 
diagram becomes more complex as the two coexistence regions partially overlap [17]. 

We consider a single bilayer membrane composed of x mole fraction of saturated lipid 
and (1 — x) of unsaturated lipid. The two lipids are taken to have different gel transition 
temperatures originating from different chain length, degree of saturation, or hydrophilic head 
group. We assume the same area per molecule for both species, and ignore lipid exchange with 
the surrounding solvent. The total free energy /^^ = ff^ + comprises (i) the free energy 
ff^ of mixing and binary interactions, and (ii) the "stretching" energy /j^, which controls 
changes in bilayer thickness. 

The free energy per site is the sum of the entropy of mixing and enthalpy. It can be 
written within a Bragg- Williams (mean-field) theory as 

/f (x) = kBT [xlogx+ (1 - x) log(l - x)] - \Jx^, (1) 

where fcs is the Boltzmann constant, T is the temperature, J > is an attractive interaction 
parameter that enhances demixing. Linear terms in x can be disregarded because they merely 
shift the origin of the chemical potential. 

To describe the gel transition, which involves chain ordering and stiffening, we introduce 
a rescaled membrane thickness ^ = {5 — Sq)/So as an order parameter, where S is the ac- 
tual membrane thickness and So is the constant membrane thickness in the disordered phase 
corresponding to the liquid phase [18]. Note that summarizes changes in various degrees 
of freedom, including the conformations of the hydrocarbon chains and their inter-chain cor- 
relations. Since the gel transition is first order, an appropriate phenomenological Landau 
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Fig. 1 - Calculated phase diagrams of binary mixtures of saturated and unsaturated lipids as a 
function of mole fraction of the saturated lipid x and temperature T for (a) J — 4.0fcBrj*, and (b) 
J = 5.0 ksT*. The other parameters are a'2 = 174fcB,a3 = -307fcBr3*,a4 = 613 ksT*, T* = 240K, 
Tg* = 260K as estimated in [18]. The gel transition temperature is Tg = 291K and 311K for x — 
and 1, respectively. The ordered (large 1/)) and disordered (1/) — 0) phases are respectively denoted 
by O and D. The critical point is indicated by a filled circle, and the triple point by Tr. The critical 
points are located at (a) Xc — 0.498, Tc = 280K, and (b) Xc = 0.5, Tc = 326K, respectively. 

expansion of the "stretching" free energy per site in powers of is [18] 

/f (x, V') = ^a',[T - T*(x)]V2 + iag^" + \a^^\ (2) 

where 02 > 0, 03 < 0, and 04 > 0. For a single component membrane (i.e., x — or 
1), T* is a reference temperature, and the first-order gel transition temperature is given 
by Tg — T* + 2a|/(9a2a4). For a binary mixture, the transition temperature depends on 
the composition x; for convenience we describe the reference temperature T*{x) as a linear 
interpolation between the two pure limits: 

T*{x)=xT: + {1-x)T:, (3) 

where T* and T* are the reference temperatures of the pure saturated and unsaturated lipids, 
respectively. With all else the same (head group size, interactions, and chain length, etc.), 
we expect T* > T*, because unsaturated lipids break up the crystallizing tendencies. Note 
that eq. Q leads to a coupling term xip'^ [19]. We have neglected a lower order bilinear term 
xip, which induces a small temperature and composition dependence in Sq. Another possible 
coupling term x'^ip simply renormalizes the interaction parameter J. 

Combining eqs. (|1I3|I . we obtain the total free energy /^^. After minimizing with respect 
to ipj the two-phase region is obtained by the Maxwell construction. In fig. ^ we show two 
typical phase diagrams of binary lipid mixtures. For small J (fig. H^a)) which corresponds 
experimentally to mixtures with weak segregation tendency, a cigar-like coexistence region 
is obtained between a disordered (D) phase where ip = (liquid) and an ordered (O) phase 
where ip > (gel). At temperatures below the cigar-shape region, there is another coexistence 
region between two ordered phases, 01-1-02, with -0i < tp2- This type of phase diagram was 
experimentally observed for DEPC-DPPC binary lipid systems [20]. When J is increased, the 
01-1-02 coexistence region extends into and beyond the D-|-0 coexistence region (fig-db)). 
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Fig. 2 - (a) Schematic phase prism of a ternary system consisting of saturated lipid (S), unsaturated 
lipid (U), and cholesterol (C). The gray and black regions on the constant temperature triangle plane 
are two-phase and three-phase coexisting regions, respectively, (b) Lattice model for a lipid-cholesterol 
mixture. Each cholesterol molecule (•) forms a dimer with a neighboring lipid molecule (o). 



In this more complex phase diagram, there are two triple points (denoted as Tr) at which 
the two disordered phases and ordered phase (D1-I-D2+02), or the disordered phase and two 
ordered phases (Dl+01+02) coexist. Many features of this phase diagram can be seen for 
DEPC-DPPE lipid mixtures which have a stronger segregation tendency because both the 
head and the tail moieties are different [20] . 

Model for lipid-cholesterol systems. - Next we discuss the role of sterols, such as choles- 
terol and lanosterol, on the phase behavior. On one hand, a small amount of cholesterol desta- 
bilizes the gel phase in favor of a liquid- disordered (Ld) phase [14, 15]. Substantial cholesterol, 
on the other hand, stabilizes a liquid- ordered (Lo) phase in which the lipid hydrocarbon tails 
are extended, but maintain high lateral mobility. This reflects the dual molecular mechanism 
of a cholesterol molecule: it can act as (i) an "impurity" and weakens the inter-lipid interac- 
tions for ordering or (ii) a "chain rigidifier" and induces conformational order in neighboring 
lipid chains [13]. Moreover, recent experiments using atomic force microscopy showed that 
the cholesterol-rich domains are thicker than the cholesterol-poor regions [21,22], supporting 
our assumption that the local membrane thickness can be taken as the order parameter. 

Based on these observations, we propose a model for lipid-cholesterol binary mixtures. 
Consider a membrane with cholesterol mole fraction c and lipid mole fraction (1 — c). There 
are three contributions to the total free energy /^'^ = f['^ + /l'^ + /|'^: (i) the entropy of mixing 
fi^{c) between lipid and cholesterol; (ii) the "stretching" energy /l'^(V') of lipid molecules; and 
(iii) direct coupling terms between lipid and cholesterol, /|'^(c, V'), which takes into account 
the effect of cholesterol on the bilayer thickness. 

We note that the coexistence region in the lipid-cholesterol system is typically found only 
for c < 0.5 [8,14,15]. For cholesterol concentration above roughly 0.5, the entire bilayer 
becomes unstable. To account for such observations we consider "condensed complexes" , fol- 
lowing [23] , and assume that each cholesterol molecule in the membrane irreversibly dimerizes 
with a single lipid molecule. Since all cholesterol molecules form dimers, as shown in fig.E^b), 
the entropy of mixing is taken between lipid-cholesterol dimers and the remaining monomeric 
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Fig. 3 - Phase diagrams of binary lipid-sterol mixtures as a function of mole fraction c of sterol 
and temperature T for (a) cholesterol, Fi — llfcBT*, F2 — 47 ksT*; and (b) lanosterol, Fi — 
28 ksT* ,T2 = 20 ksT*. The other parameters are the same as in fig. H except T* = 260K. The gel 
transition temperature at c = is Tg = 311K. The ordered (large i/j) and disordered ( non-zero small 
1/)) phases are respectively denoted by O and D. The critical point indicated by a filled circle in (a) is 
located at Cc = 0.492, Tc = 338K. 



lipids, and is given by a Flory-Huggins form: 

/f=(c) = kBT[c\og2c+il - 2c)log(l - 2c)] . (4) 

This free energy is valid only for < c < 0.5, since each cholesterol molecule is paired with a 
neighboring lipid. The first term is the entropy of dimers having an area fraction of 2c, while 
the second term accounts for the entropy of lipid monomers of area fraction (1 — 2c) . Although 
extensions of the model to include excess free cholesterol and/or formation of trimers, etc. is 
possible [24], the present situation is enough to describe the qualitative phase behavior. 

The Landau free energy describing the structural phase transition of the membrane is 
given by /l'^(V') and is the same as /|^ of eq. 

f2{^) = h'2{T - r*)^' + iflsV'' + ia4V^ (5) 

where the order parameter ip is again the relative bilayer thickness, and T* is the reference 
temperature of the pure lipid system. 

The simplest free energy to account for the previously mentioned dual effects of cholesterol 
can be expressed by the following phenomcnological coupling terms: 

/3^^(c,v) = iricV'-ir2cV, (6) 

where Fi > and r2 > are the coupling constants. The first term expresses the fact that 
a small amount of cholesterol (c > 0) acts as an "impurity" . It interferes with the crystalline 
ordering, and enhances disordered chains (smaller ■(/;). This term is necessary because, as 
shown in the experimental phase diagram [14,15], the gel-L^ coexistence temperature decreases 
upon adding cholesterol to the pure lipid system. The second term corresponds to the "chain 
rigidifying" effect reflecting the fact that a larger amount of cholesterol favors ordered tail 
states and hence larger Meanwhile, this coupling induces lipid-cholesterol phase separation 
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since it is proportional to with a negative coefScient. In other words, the effective hpid- 
cholesterol interaction depends on the conformational states of the lipid chains. 

Adding eqs. H4I6|I . we first minimize /^'^ with respect to i/) and then construct the phase 
diagrams as shown in fig. |31 In fig. Ola), there are three different coexisting regions. Due to 
the Fi coupling term in eq. ©, the region of the disordered (D) phase, characterized by a 
non-zero small value of i/;, widens at expense of ordered (01) phase, characterized by a large 
value of -0, upon the addition of a small amount of cholesterol (c < 0.1). Note that D- and 01- 
phases respectively correspond to Ld and gel phases in the experiments. For larger c, however, 
the F2 coupling term overcomes the first term, and the region of the D-phase narrows in favor 
of 02-phase (c > 0.1), which corresponds to the Lo phase in which the order parameter ^jj 
takes the value much larger than that in the D-phase. The obtained phase diagram for lipid- 
cholesterol systems agrees semiquantitatively with the experimental one [14,15]. The D- and 
02-phases have the same symmetry and are continuously connected, with a critical point at 
Cc = 0.492 < 0.5. The appearance of the critical point can be understood by noticing that 
the coupling terms in eq. ((BJ are linear in ip, and act as an external field coupled io ip. In 
general, the first-order transition becomes continuous when the applied external field is strong 
enough [25]. 

For larger Fi and smaller r2 we obtain the phase diagram shown in fig. Ofb). In this 
case, the effective external field is too weak to eliminate the first-order transition, and the 
critical point does not exist. This phase diagram resembles that of binary lipid-lanosterol 
mixtures [26]. Compared to cholesterol, lanosterol has three additional methyl groups and 
has a structurally rougher hydrophobic part, leading to a weaker enhancement of lipid tail 
stretching (smaller F2). Moreover, an extra double bond in the lanosterol hydrocarbon tail 
can be expected to inhibit crystallinity (larger Fi). In fact, lanosterol is considered to be a 
precursor to cholesterol in the evolutionary pathway [26]. Although cholesterol can stabilize 
the region of coexistence between the D- and 02-phases [i.e., Ld+Lo) as shown in fig. Efa), 
such a coexistence has not been identified for the lipid-lanosterol mixtures [26] . The fact that 
the coexistence region ends at c = 0.5 is due only to the artificial truncation of the model 
at c — 0.5. If free excess lanosterol were allowed, the phase boundaries would continue until 
the membrane lost its integrity. However, such a break up of the membrane is not considered 
here. 

Discussion. ~ Within a few simplifying assumptions we have obtained phase diagrams 
that successfully mimic the experimental ones. A few points merit further discussion. We used 
only one order parameter ip (related to local membrane thickness) to describe the membrane 
internal degrees of freedom of chain ordering and crystallinity. This assumption is well moti- 
vated for pure lipid membrane where both ordering and crystallization occur simultaneously. 
Although addition of cholesterol may separate the chain ordering and freezing transitions and 
stabilize the Lo phase, the thickness ip is sufficient here to distinguish between the gel (01), 
Lo (02), and Ld (D) phases. Note that the symmetry of the D- and 02-phases is the same 
(hence the critical point), and differs from the solid- like and local crystalline order of the 
01-phase. As an extension to the present model, it may be of interest to use two distinct 
order parameters. 

Finally, we discuss the qualitative behavior of the full ternary phase diagram for saturated 
(S) and unsaturated (U) lipids in the presence of cholesterol (C). The phase diagrams in 
figs.^a) and|2Ia) correspond to the two side views of the schematic phase prism in fig.|2a), 
namely the US-T and the SC-T planes, respectively. A similar phase diagram to fig.EIa) is ex- 
pected for the UC-T plane. Consider a cut of the phase prism at a fixed temperature as shown 
in fig.|2Ia). Both the D-l-0 coexistence region on the US-T plane and the 01-1-02 coexistence 
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region on the SC-T plane extend onto this triangle plane, and may form a three-phase coexist- 
ing region (D-|-01-|-02, i.e., Ld+gel+Lo) when they meet each other. If the phase separation 
does not occur between U and C for this temperature, a critical point between D and 02 is 
expected to appear. As the temperature is increased, the line of critical points may end when 
the three-phase region disappears. These arguments are consistent with recent experimental 
studies of the phase morphology inside the phase prism [6-9] . Further calculations of the full 
ternary phase diagrams are in progress. 

In summary, we have analyzed theoretically saturated and unsaturated lipid-lipid mixtures, 
and lipid- cholesterol mixtures. The coupling between composition and the internal degree of 
freedom (structural chain melting transition) allows us to obtain phase diagrams that agree 
well with experiment, and to explain some of the major features leading to creation of rafts. 
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